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1. Introduction 

A prominent theme of modern condensed matter phy- 
sics is electronic transport - in particular the electrical 
conductivity - of disordered metallic systems at very 
low temperatures. From the Landau theory of weakly 
interacting Fermi liquids one expects the essential as- 
pects of the situation to be captured by the single-elec- 
tron approximation. Mathematical models that have 
been proposed and studied in this context include ran- 
dom Schrodinger operators and band random matrices. 

If the physical system has infinite size, two distinct 
possibilities exist: the quantum single-electron motion 
may either be bounded or unbounded. In the former 
case the disordered electron system is an insulator, in 
the latter case a metal with finite conductivity (if the 
electron motion is not critical but diffusive). Metallic 
behavior is expected for weakly disordered systems in 
three dimensions; insulating behavior sets in when the 
disorder is increased or the space dimension reduced. 

The main theoretical tool used in the physics liter- 
ature on the subject is the supersymmetry method pi- 
oneered by Wegner and Efetov (1979-1983). Over the 
past twenty years, physicists have applied the method in 
many instances, and a rather complete picture of weakly 
disordered metals has emerged. Several excellent re- 
views of these developments are available in print. 

From the perspective of mathematics, however, the 
method has not always been described correctly, and 
what is sorely lacking at present is an exposition of how 
to implement the method rigorously. [Unfortunately, 
the correct exposition by Schafer and Wegner (1980) 
was largely ignored or forgotten by later authors.] In 
this encyclopedia article an attempt will be made to help 
remedy the situation, by giving a careful review of the 
Wegner-Efetov supersymmetry method for the case of 
Hermitian band random matrices. 



Let V be a unitary vector space of finite dimension. 
A Hermitian random matrix model on V is defined by 
some probability distribution on Herm(V), the Hermi- 
tian linear operators on V. You may fix some orthonor- 
mal basis of V and represent the elements H of Herm(V ) 
by Hermitian square matrices. 

Quite generally, probability distributions are charac- 
terized by their Fourier transform or characteristic func- 
tion. In the present case this is 

where the Fourier variable K is some other linear op- 
erator on V, and (...) denotes the expectation value 
w.r.t. the probability distribution for H. Later it will 
be important that, if £l(K) is an analytic function of K, 
the matrix entries of K need not be from R or C but can 
be taken from the even part of some exterior algebra. 

The probability distributions to be considered in this 
article are Gaussian with zero mean, (H) = 0. Their 
Fourier transform is also Gaussian: 

with J some quadratic form. We now describe / for 
a large family of hierarchical models that includes the 
case of band random matrices. 

Let V be given a decomposition by orthogonal vector 
spaces: 

V = Vi@V 2 @ ...@V\ A \ . 

You should imagine that every vector space V, corre- 
sponds to one site i of some lattice A, and the total num- 
ber of sites is |A|. For simplicity, we take all dimen- 
sions to be equal: dimVi = . . . = dim^ A | = N. Thus 
the dimension of V is N\A\. The integer N is called the 
number of orbitals per site. 

If n, is the orthogonal projector on the linear sub- 
space Vj CV, we take the bilinear form / to be 

|A| 

J{K,K') = £ JiMUiKUjK') , 

where the coefficients Jjj are real, symmetric, and pos- 
itive. This choice of J implies invariance under the 
group of unitary transformations in each subspace: 

q/ =u(Vi) xu(v 2 ) x ■••xU(Vj A |) • 

Clearly, D.(K) = Q.(UKU^ 1 ) or, equivalently, the prob- 
ability distribution for H is invariant under conjugation 
H i-> UHU~ X , for U € . 

If {ef } a =i,...,jv is an orthonormal basis of V,-, we de- 



fine linear operators E"^ : Vj — > V, by E\ 



jahJb 



ey. 



By 



evaluating J(Eff t E$f) = J ij 8 U! 8 jf 8 aa! 8 W one sees that 
the matrix entries of H are all statistically independent. 
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By varying the lattice A, the number of orbitals N, 
and the variances Jjj, one obtains a large class of Hermi- 
tian random matrix models, two prominent subclasses 
of which are the following: 

1. For |A| = 1, one gets the Gaussian Unitary En- 
semble (GUE). Its symmetry group is ^ = U(iV), 
the largest one possible in dimension N = dimV. 

2. If \i — j\ denotes a distance function for A, and / 
a rapidly decreasing positive function on M + of 
width W, the choice J u = f(\i- j\) with N = 1 
gives an ensemble of band random matrices with 
band width W and symmetry group % = U(1)I A I . 

Beyond being real, symmetric and positive, the vari- 
ances Jij are required to have two extra properties in 
order for all of the following treatment to go through: 

• They must be positive as a quadratic form. This 
is to guarantee the existence of an inverse, which 
we denote by wtj = (J~ l )ij. 

• The off-diagonal matrix entries of the inverse must 
be non-positive: Wjj < for i / j. 



3. Basic tools 

3.1. Green's functions. — A major goal of random- 
matrix theory is to understand the statistical behavior 
of the spectrum and the eigenstates of a random Hamil- 
tonian H. Spectral and eigenstate information can be 
extracted from the Green's function, i.e. from matrix 
elements of the operator (z — H)~ l with complex pa- 
rameter z 6 C\K. For the models at hand, the good 
objects to consider are averages of "1/ -invariant observ- 
ables such as 

(1) G; i) (z) = <Trn i -(z-//)- 1 >, 

(2) G^(z u zi) = (Trn^zi -H)- l TVj{z 2 -H)- 1 ) . 

The discontinuity of G^p (z) across the real z-axis yields 
the local density of states. In the limit of infinite vol- 

(2) 

ume (|A| —> oo), the function G)j (zi,zi) for z\ =E+ie, 
Z2=E — ie, real energy E, and e > going zero, gives 
information on transport, e.g. the electrical conductivity 
by the Kubo-Greenwood formula. 

(2) 

Mathematically speaking, if G)j(E + ie,E — ie) is 
bounded (for infinite volume) in e and decays alge- 
braically with distance | i — j\ at e = 0+, the spectrum is 
absolutely continuous and the eigenstates are extended 
at energy E. On the other hand, a pure point spectrum 
and localized eigenstates are signalled by the behavior 

G^ ~ £~ l e~^'~i\ with positive Lyapunov exponent A. 



3.2. Green's functions from determinants. — For any 

pair of linear operators A, B on a finite-dimensional vec- 
tor space V, the following formula from basic linear al- 
gebra holds if A has an inverse: 



— Det(A + tB) 
at 



(=0 



Det(A)Tr(A _1 B) 



Using it with A = z — H and z G C \ R, all Green's func- 
tions can be expressed in terms of determinants; for ex- 
ample, G\f(w,z) = 



N 
I 



d 2 /Det(w-H)Det(z-H + tE. 



s=t=0 



J^ =x dsdt \ Det(w -H- sEf)T>Qi(z - H) 

It is clear that, given a formula of this kind, what one 
wants is a method to handle ensemble averages of ratios 
of determinants. This is what's reviewed in the sequel. 

3.3. Determinants as Gaussian integrals. — Let the 

Hermitian scalar product of the unitary vector space V 
be written as <pi, <j?2 i— > ($i, (fa), and denote the adjoint 
or Hermitian conjugate of a linear operator A on V by 
A*. If 9\tA := i(A +A*) > 0, the standard Lebesgue 
integral of the Gaussian function <p ^(vAv) makes 
sense and gives 



(3) 



j e -(«u«>) =D etA- 1 , 



where it is understood that we are integrating with the 
Lebesgue measure on (the normed vector space) V nor- 
malized by / e~^ i<? ) = 1. The same integral with anti- 
commuting y instead of the (commuting) (p G V gives 



(4) 



-(<MV) = DetA . 



This basic formula from the field theory of fermionic 
particles is a consequence of the integration over anti- 
commuting variables actually being differentiation: 

d 2 

dijtidYiftyuYi,...) := f(%,y u ...) . 



/■ 



4. Fermionic variant 



The supersymmetry method of random-matrix the- 
ory is a theme with many variations. The first variation 
to be described is the "fermionic" one. To optimize the 
notation, we now write djXNj{H) for the density of the 
Gaussian probability distribution of H: 



(F(H)}= J F(H)dli N j{H) . 

All determinants and traces appearing below will be 
taken over vector spaces that are clear from the context. 



THE SUPERSYMMETRY METHOD OF RANDOM MATRIX THEORY 



3 



Let zi, ■ ■ ■ ,z n be any set of n complex numbers, put 
z := diag(zi , . . . ,z n ) for later purposes, and consider 

(5) OffifaJ) = ffl Det( Za - H) dn N ,j(H) . 

■' ce=l 

The supersymmetry method expresses this average of a 
product of determinants in an alternative way, by inte- 
grating over a "dual" measure as follows. 

Introducing an auxiliary unitary vector space C", one 
associates with every site i of the lattice A an object 
Qi £ Herm(C"), the space of Hermitian n x n matri- 
ces. If dQi for i = 1, . . . , |A| are Lebesgue measures on 
Herm(C"), one puts DQ = const x YljdQj and 

(6) dv n j(Q) := e-5^( r '^ Tre %2 . 

The multiplicative constant in DQ is fixed by requiring 
the density to be normalized: / dv n j{Q) = 1. By com- 
pleting the square, this Gaussian probability measure 
has the characteristic function 

f eVWi dv n j(Q) =e-^A™' , 

where the Fourier variables K\,. .. ,K\&\ are n x n matri- 
ces with matrix entries taken from C or another com- 
mutative algebra. 

The key relation of the fermionic variant of the su- 
persymmetry method is that the expectation of the prod- 
uct of determinants © has another expression as 




(i = v— T). The strategy of the proof is quite sim- 
ple: one writes the determinants in both expressions for 
£ljj e ™ as Gaussian integrals over nAf |A| complex fermio- 
nic variables Yi > • • • > Yn (each Ya is a vector in V with 
anti-commuting coefficients), using the basic formula 
©. The integrals then encountered are essentially the 
Fourier transforms of the distributions dflNj(H) resp. 
dv n j(Q). The result is 

J e - Ey Zr( Vr, Vr) e - \ Ey hj lap (Va A Vp) ifp Va) 

for both expressions of ^Jj e ™. In other words, although 
the probability distributions dpL^j{H) and dv n j{Q) are 
distinct (they are defined on different spaces), their char- 
acteristic functions coincide when evaluated on the Fou- 
rier variables K = Y,a¥a(Ya,') for H and (Kj) a p = 
(yir a ,rLy/fi) for Qj. This establishes the claimed equal- 
ity of the expressions © and © for Q.^ r ™(z,J). 

What's the advantage of passing to the alternative 
expression by dv n j(Q)7 The answer is that, while H 
is made up of independent random variables, the new 
variables Qi, called the Hubbard-Stratonovich field, are 



correlated: they interact through the "exchange" con- 
stants Wjj = (J^ l )ij. If that interaction creates enough 
collectivity, a kind of mean-field behavior results. 

For the simple case of GUE (|A| = 1, w u =N/X 2 ) 
with z\ = ■ ■■ = z n =E, one gets the relation 

(Det" (E-H)) = J Det N (E-iQ)e~W TT &dQ , 

the right-hand side of which is easily analyzed by the 
steepest descent method in the limit of large N. 

For band random matrices in the so-called ergodic 
regime the physical behavior turns out to be governed 
by the constant mode Q\ = ■■■ = Qui - a fact that can 
be used to establish GUE universality in that regime. 

5. Bosonic variant 

The bosonic variant of the present method, due to 
Wegner, computes averages of products of determinants 
placed in the denominator: 

(8) a b n ^(z,J) = ff\ Det-\z a -H)diiNj(H) , 

J a=\ 

where we now require 3m Za ^ for all a = 1,. . . ,n. 
Complications relative to the fermionic case arise from 
the fact that the integrand in (J8jl has poles. If one re- 
places the anti-commuting vectors y a by commuting 
ones q> a , and then simply repeats the previous calcula- 
tion in a naive manner, one arrives at 




where the integral is still over Qj G Herm(C"). The 
calculation is correct, and relation © therefore holds 
true, provided that the parameters z\ , . ■ • ,Z n all lie in the 
same half (upper or lower) of the complex plane. To ob- 
tain information on transport properties, however, one 
needs parameters in both the upper and lower halves; 
see the paragraph following ©. The general case to 
be addressed below is 3mz a > for a = 1, . .. ,p, and 
3mz a < for a = p + 1, . .. ,n. Careful inspection of 
the steps leading to equation © reveals a convergence 
problem for < p < n. In fact, © with Qj in Herm(C) 
turns out to be false in that range. Learning how to re- 
solve this problem is the main step toward mathematical 
mastery of the method. Let us therefore give the details. 

If s a := sgnCJmza, the good (meaning convergent) 
Gaussian integral to consider is 

j ^ a s a [9a,(i a -H)<p a ) = flDet- l (-is a (z a -H)) . 
J a=\ 

To avoid carrying around trivial constants, we now as- 
sume i(" _2 P) A 'l A l = l. Use of the characteristic function 
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of the distribution for H then gives 

(10) QE&M = j e^rW^-Pr) 

The difficulty of analyzing this expression stems from 
the "hyperbolic" nature (due to the indefiniteness of the 
signs s a = ±1) of the term quartic in the (p a ,(f>a- 

5.1. Fyodorov's method. — The integrand for £2 bos is 
naturally expressed in terms of n x n matrices M, with 
matrix elements (Mi) a p = (<p a ,rL<pft). These matrices 
lie in Herm + (C"), i.e. they are non-negative as well as 
Hermitian. Fyodorov's idea was to introduce them as 
the new variables of integration. To do that step recall 
the basic fact that, given two differentiable spaces X 
and Y and a smooth map y : X — ► Y, a distribution jj. 
on X is pushed forward to a distribution y/"(jLt) on 7 by 
y(ju) [/] := M t/° V]> where / is any test function on 7. 
We apply this universal principle to the case at hand 

by identifying X with V n , and Y with (Herm+(C")) |A| , 
and \j/ with the mapping that sends 

(<p b ...,<p n ) eX to (Mi,...,M| A |) £Y 

by (Afi)«p = (9a,n,-9p)- On X = V" we are integrat- 
ing with the product Lebesgue measure normalized by 
J Q-L a {9a,9a) = i We now want the push forward of 
this flat measure (or distribution) by the mapping \y. In 
general, the push forward of a measure is not guaran- 
teed to have a density but may be singular (like a Dirac 
8 -distribution). This is in fact what happens if N < n. 
The matrices M, then have less than the maximal rank, 
so they fail to be positive but possess zero eigenvalues, 
which implies that the flat measure on X is pushed for- 
ward by y mto the boundary of Y. For N >n, on the 
other hand, the push forward measure does have a den- 
sity on Y; and that density is ]jj=i (petMi) N ~ n dMu as 
is seen by transforming to the eigenvalue representation 
and comparing Jacobians. The dM, are Lebesgue mea- 
sures on Herm(C"), normalized by the condition 

/ tr TtM '(petMif- n dM i = [ e -E«(fc" n '*») = 1 . 

JMj>0 J 

Assembling the sign information for 3mz, a in a diag- 
onal matrix s := diag(si, . . . ,s n ), and pushing the inte- 
gral over X forward to an integral over Y with measure 
DM := YlidMi, we obtain Fyodorov's formula: 

(11) n^( z ,y) = J e -hi-ijJiM^ t sM } ) 

x e L k Tr(iszM k +(N-n)lnM k ) DM _ 

This formula has a number of attractive features. One 
is ease of derivation, another is ready generalizability 



to the case of non-Gaussian distributions. The main 
disadvantage of the formula is that it does not apply 
to the case of band random matrices (because of the re- 
striction N > n); nor does it combine nicely with the 
fermionic formula Q to give a supersymmetric formal- 
ism, as one formula is built on Jjj and the other on Wij. 

Note that (TTTt clearly displays the dependence on 
the signature of 3m z '■ you cannot remove the si,...,s„ 
from the integrand without changing the domain of in- 
tegration Y = (Herm + (C")) ' A '. This important feature 
is missing from the naive formula 

Setting q = n — p, let \J(p,q) be the pseudo-unitary 
group of complex nxn matrices T with inverse = 
sT*s. Since |DetT| = 1 for T G U(p,q), the integration 
domain Y and density DM = Y[jdMj of Fyodorov's for- 
mula are invariant under U(p,q) transformations M; i— > 
TMjT*, and so is actually the integrand in the limit 
where all parameters z\,---,z n become equal. Thus the 
elements of \J(p,q) are global symmetries in that limit. 
This observation holds the key to another method of 
transforming the expression (flOl . 

5.2. The method of Schafer and Wegner. — To res- 
cue the naive formula (O, what needs to be abandoned 
is the integration domain Herm(C") for the matrices Qj. 
The good domain to use was constructed by Schafer and 
Wegner, but was largely forgotten in later physics work. 

Writing {Mu) a p = (q> a ,^-k(pp) as before, consider 
the function 

(12) F M (Q) = e^ i J w, j Tt ( sQi+ ^( sQ j + ^~'^ TTMkQk , 
viewed as a holomorphic function of 

G = (Gi,.--,e|A|)GEnd(C") |A| • 

If the Gaussian integral / Fm(Q)DQ with holomorphic 
density DQ = fT ; dQi is formally carried out by com- 
pleting the square, one gets the integrand of (110b . This 
is just what we want, as it would allow us to pass to a 
<2-matrix formulation akin to the one of Section 0] But 
how can that formal step be made rigorous? To that end, 
one needs to (i) construct a domain on which |^m(<2)| 
decays rapidly so that the integral exists, and (ii) justify 
completion of the square and shifting of variables. 

To begin, take the absolute value of Fm{Q)- Putting 
\{Qj + Q)) =■ KeQj and ^{Q } - Q*) =: 3m Qj, you 
have \F M \ = erk(h+h+h) with 

fi (fi) = £ wyTr (s 3m Q x + z) (s 3m Qj + z)+c.c, 

U 

MQ) = -2£wyTr(j3teGf)(jJReGy) , 

ij 

MQ) =4£Tr (Mi + slmzEjWij) VKtQi . 
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These expressions suggest making the following choice 
of integration domain for 2, (/ = 1, . . . , |A|). Pick some 
real constant A > and put 

f Ff 
P 



with T t G V(p,q), P+ G Herm(O'), Pr G Herm(C«). 
The set of matrices G; so defined is referred to as the 
Schafer-Wegner domain X?' q . The range of the field 

Q = (Gl , • • • , Q\A\) is the direct product X := (xg 9 ) |A| . 

To show that this is a good choice of domain, we first 
show convergence of the integral J^- F M (Q)DQ. First 
of all, the matrices P; commute with s, so 



/i(G) 



x 



2 9ie £ w u Tr (P* + sz) (Pj + sz) 



Since the coefficients Wij are positive as a quadratic 
form, this expression is convex (with a positive Hes- 
sian) in the Hermitian matrices p. Second, the function 

-2A 2 £w i7 Tr (P^T 1 



/ 2 (G) 



is bounded from below by the constant — 2A 2 w£ ( -w,7. 
This holds true because w,y is negative for i ^ j, and 
because T t T* > and the trace of a product of two pos- 
itive Hermitian matrices is always positive. Third, 

MQ) = 4A £Tr (Mi + sJvaz I y w, 7 ) T(T* 

is positive, as (...) is positive Hermitian. As long as 
s 3mz > the function f$ goes to infinity for all possible 
directions of taking the 7} to infinity on\J(p,q). 

Thus when the matrices G; are taken to vary on the 
Schafer-Wegner domain X^ q , the absolute value \Fm\ = 
Q-iifi+fz+A) decays rapidly at infinity. This establishes 
the convergence of f%- Fm(Q)DQ. 

Next, let us count dimensions. The mapping T i— > 
PP* for T G U(p,q) =: G is invariant under right mul- 
tiplication of T by elements of the unitary subgroup 
H := U(/j) x \J(q) - it is called the Cartan embedding 
of G/H into G. The real manifold G/H has dimension 
2/?g and so does its image under the Cartan embedding. 
Augmenting this by the dimension of Herm(C p ) and 
Herm(C^) (from p), one gets dimX^' q = 2pq + p 2 + 
q 2 = (p + q) 1 = which is as it should be. 

Finally, why can one shift variables and do the Gaus- 
sian integral over Q (with translation-invariant DQ) by 
completing the square? This question is legitimate as 
the Schafer-Wegner domain X^' q lacks invariance under 
the required shift, which is G; l— ► Qi — isz + E; /y *M/ J. 

To complete the square in (fl2l . introduce a parame- 
ter ? G [0, 1 [ and consider the family of shifts 

Qi ^Qi + t(-isz + Zj Jij sMjs) . 



For fixed t, this shift takes X = (X?' q ) into another 
domain, 3£(t). Inspection shows that the function (fl2l 
still decays rapidly (uniformly in the M/) on 3£(f), as 
long as t < 1 . Without changing the integral one can add 
pieces to 3£(t) (for t < 1) at infinity to arrange for the 
chain X — 3C(f) to be a cycle. Because 5£{t) is homo- 
topic to JT(0) = 5£ , this cycle is a boundary: there ex- 
ists a manifold of dimension dim X + 1 such that 
dW(t) = X - %{t). Viewed as a holomorphic dif- 
ferential form of degree (« 2 |A|,0) in the complex space 
End(C")l A l , the integrand CO := F M (Q)DQ is closed (i.e. 
dco = 0). Therefore, by Stokes' theorem, 

co- co= co= d&) = 0, 

JX JX(t) JdW{t) J»{t) 

which proves J S -^F M (Q)DQ = $ %■ F M (Q)DQ, inde- 
pendent of t. (This argument does not go through for 
the non-rigorous choice sQi := TiPjTf usually made!) 
In the limit t — ► 1 , one encounters the expression 



/ F M (Q)DQ= f dv n j{iQ) 

JX(l) JX 

x g— I Ey Ji/Tr (sMiSMj) Tr (szA4) 



with <iv„ j as in ©. The normalization integral over 
<5T is defined by taking the Hermitian matrices p to be 
the inner variables of integration. The outer integrals 
over the 7} then demonstrably exist, and one can fix 
the (otherwise arbitrary) normalization of DQ by setting 
/ % dv n j (iQ) = 1 . Making that choice, and comparing 
with <fT0b . one has proved 



/-vbos 



X 



F (M,) a p = (faM,<p p )(Q) D Q 



The final step is to change the order of integration 
over the Q- and (^-variables, which is permitted since 
the G-integral converges uniformly in (p. Doing the 
Gaussian (^-integral and shifting Qk — > Qk~ isz, one 
arrives at the Schafer-Wegner formula for ^ c 



bos . 



(13) Q^^w- 1 ) = J &U"i}*{!>Q0Q) 



x e 



-JVLtTrln(e t -i5 Z ) 



DQ, 



which is a rigorous version of the naive formula ®. 
Compared to Fyodorov's formula, it has the disadvan- 
tage of not being manifestly invariant under global hy- 
perbolic transformations Gi l— ► TQjT* (the integration 
domain 3C isn't invariant). Its best feature is that it 
does apply to the case of band random matrices with 
one orbital per site (N = 1). 
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6. Supersymmetric variant 

We are now in a position to tackle the problem of av- 
eraging ratios of determinants. For concreteness, we 
shall discuss the case where the number of determi- 
nants is two for both the numerator and the denomi- 
nator, which is what is needed for the calculation of the 
function (zi,Z2) defined in equation We will 
consider the case of relevance for the electrical conduc- 
tivity: z\ = E + i£, Z2 = E - ie, with E G R and e > 0. 

(21 

A <2-integral formula for - (zi,Z2) can be derived 
by combining the fermionic method for 

(Det( Zl - H) Det(z 2 -H + t 2 E?f)) 

with the S chafer- Wegner bosonic formalism for 

(Der \ Zl -H-hEf) Der 1 (z 2 - H) ) , 

and eventually differentiating with respect to t\,t% at 
t\ = ?2 = and summing over a,b; see Section 13721 All 
steps are formally the same as before, but with traces 
and determinants replaced by their supersymmetric ana- 
logs. Having given a great many technical details in 
Sections |4]and|5j we now just present the final formula 
along with the necessary definitions and some indica- 
tion of what are the new elements involved in the proof. 

Let <2bb, Qff, Qbf and <2fb each stand for a 2 x 2 
matrix. If the first two matrices have commuting entries 
and the last two anti-commuting ones, they combine to 
a 4 x 4 superniatrix: 

( <2bb <2bf 
\<2fb 2ff, 

Relevant operations on supermatrices are the supertrace, 

STr£ = TreBB-Tre FF , 
and the superdeterminant, 

SDet 6 = - D ^ l . 
Det(<2pp — <2fb2bb 2bf) 

These are related by the identity SDet = exp oSTr o In 
whenever the superdeterminant exists and is nonzero. 

In the process of applying the method described ear- 
lier, a supermatrix Qj gets introduced at every site i of 
the lattice A. The domain of integration for each of 
the matrix blocks (2,)bb Q = 1 5 • • • , |A|) is taken to be 
the S chafer- Wegner domain X.' (with some choice of 
A > 0); the integration domain for each of the (<2,-)ff is 
the space of Hermitian 2x2 matrices, as before. 

Let £gg be the 4x4 (super)matrix with unit entry 
in the upper left corner and zeroes elsewhere; simi- 
larly, £pp has unity in the lower right corner and ze- 
roes elsewhere. Putting s = diag ( 1 , — 1,1,1) and z = 
diag(zi ,Z2,Z\ ,Zt), the supersymmetric <2-integral for- 

(2) 

mula for the generating function of G\- - obtained by 



combining the Schafer- Wegner bosonic method with the 
fermionic variant - is written as 



Q 



(14) 



Det(zi - H) Det(z 2 -H + t 2 E?f) 
^ Det(zi -H- hEf) Det(z 2 - H) 



OQQl^kl W« STr (sQ k sQi) 



xe 



where the second supertrace includes a sum over sites 
and orbitals, and on setting t\ = t% = becomes 

^STrMc-te) =T isDet- N (Qr-isz) ■ 

The superintegral 'measure' DQ = Y\ r DQ r is the flat 
Berezinform, i.e. the product of differentials for all the 
commuting matrix entries in (<2,-)bb and {Q v )yf, times 
the product of derivatives for all the anti-commuting 
matrix entries in (<2>-)bf and (<2/-)fb- 

To prove the formula dl4t . two new tools are needed, 
a brief account of which is as follows. 

6.1. Gaussian superintegrals. — There exists a su- 
persymmetric generalization of the Gaussian integra- 
tion formulas given in Section EQ1 if A,D (B,C) are 
linear operators or matrices with commuting (resp. anti- 
commuting) entries, and 9^eA > 0, one has 



SDet 



B 
D 



Verification of this formula is straightforward. Using 
it, one writes the last factor in (fl4l as a Gaussian su- 
perintegral over four vectors: (f>\, (fa, y\, and \ff2- The 
integrand then becomes Gaussian in the matrices Q r . 

6.2. Shifting variables. — The next step in the proof 
is to do the 'Gaussian' integral over the supermatrices 
Q r . By definition, in a superintegral one first carries out 
the Fermi integral, and afterwards the ordinary integra- 
tions. The Gaussian integral over the anti-commuting 
parts (<2 r )BF and (<2r)FB is readily done by complet- 
ing the square and shifting variables using the fact that 
fermionic integration is differentiation: 

j^ /( £-£')=^/(s-a= jdtm. 

Similarly, the Gaussian integral over the Hermitian ma- 
trices (<2r)pF is done by completing the square and shift- 
ing. The integral over (2 r )BB» however, is not Gaus- 
sian, as the domain is not R" but the Schafer- Wegner 
domain. Here, more advanced calculus is required: these 
integrations are done by using a supersymmetric change- 
of-variables theorem due to Berezin to make the neces- 
sary shifts by nilpotents. (There is not enough space 
to describe this here, so please consult Berezin's book.) 
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Without difficulty one finds the result to agree with the 
left-hand side of Eq. (fT4b . thereby establishing that for- 
mula. 



7. Approximations 

All manipulations so far have been exact and, in fact, 
rigorous (or can be made so with little extra effort). 
Now we turn to a sequence of approximations that have 
been used by physicists to develop a quantitative un- 
derstanding of weakly disordered quantum dots, wires, 
films etc. While physically satisfactory, not all of these 
approximations are under full mathematical control. We 
will briefly comment on their validity as we go along. 

7.1. Saddle-point manifold. — We continue to con- 

(2) 

sider G)j (E + ie , E — ie ) and focus on E = (the center 
of the energy band) for simplicity. By varying the ex- 
ponent on the right-hand side of (fPfl i at t\ = ?2 = 0, one 
gets the following equation: 

j 

which is called the saddle-point equation. 

Let us now assume translational invariance, wu = 



f(\i — Then, if A = ^JN /Y,j w ij » the saddle-point 
equation has /-independent solutions of the form 



where for <7ff there are three possibilities: two isolated 
points ^ff = ±1 (unit matrix) coexist with a manifold 

cos 0i sin 0i e* 
, sin 0! e -1 * 1 - cos 0i 



(15) 



<7ff 



(16) 



which is 2-dimensional, whereas the solution space for 
<7bb consists of a single connected 2-manifold: 

cosh 0o sinh 0o e 1 ^ 1 
v sinh0 o e _i0 ° cosh0 o 

The solutions <7ff = ±1 are usually discarded in the 
physics literature. (The argument is that they break su- 
persymmetry and therefore get suppressed by fermionic 
zero modes. For the simpler case of the one-point func- 
tion Q and in three space dimensions, such suppres- 
sion has recently been proved by Disertori, Pinson and 
Spencer.) Other solutions for ^bb are ruled out by the 
requirement £Ke Q\ > for the Schafer-Wegner domain. 

The set of matrices (fl6l and (fT3T > - the saddle-point 
manifold - is diffeomorphic to the product of a two- 
hyperboloid H 2 with a two-sphere S 2 . Moving along 
that manifold M := H 2 x S 2 leaves the <2-field integrand 
(Tl4l unchanged (for z\ = Zi = h = ti = 0). 



One can actually anticipate the existence of such a 
manifold from the symmetries at hand. These are most 
transparent in the starting point of the formalism as given 
by the characteristic function (e.~ ljTHK ) with 

K = <pi ® <pi - (p2 <8> <P2 + % ® Vi + ¥2 <8> ¥2 ■ 

The signs of this quadratic expression are what is en- 
coded in the signature matrix s = diag(l, — 1, 1,1) (re- 
call that the first two entries are forced by 3mz\ > and 
3mz2 < 0). The quadratic form K is invariant under the 
product of two Lie groups: U(l, 1) acting on the <p's, 
and U(2) acting on the i/a's. This invariance gets trans- 
ferred by the formalism to the Q-side; the saddle-point 
manifold M is in fact an orbit of the group action of 
G := U(l, 1) x U(2) on the 0-field. In the language of 
physics, the degrees of freedom of M correspond to the 
Goldstone bosons of a broken symmetry. 

K also has a number of supersymmetries, mixing <p's 
with y's. At the infinitesimal level, these combine with 
the generators of G to give a Lie superalgebra of sym- 
metries g := u(l , 1 12). One therefore expects some kind 
of saddle-point supermanifold, say on the <2-side. 

j$ can be constructed by extending the above solu- 
tion go := diag(^BBi^FF) of the dimensionless saddle- 
point equation sqs = q~ l to the full 4x4 supermatrix 



space. Putting q = qo + q\ with q\ 
linearizing in q\, one gets 





<7fb 



<7bf 



, and 



(17) 



sqi s 



The solution space of this linear equation for q\ has di- 
mension four for all qo G M. Based on it, one expects 
four Goldstone fermions to emerge along with the four 
Goldstone bosons of M. 

For the simple case under consideration, one can in- 
troduce local coordinates and push the analysis to non- 
linear order, but things get quickly out of hand (when 
done in this way) for more challenging, higher-rank 
cases. Fortunately, there exists an alternative, coordinate- 
independent approach, as the mathematical object to be 
constructed is completely determined by symmetry! 



7.2. Riemannian symmetric superspace. — The lin- 
ear equation (TlTt associates with every point x G M a 
four-dimensional vector space of solutions V x . As the 
point x moves on M the vector spaces V x turn and twist; 
thus they form what is called a vector bundle V over M. 
(The bundle at hand turns out to be non-trivial, i.e. there 
exists no global choice of coordinates for it.) 

A section of V is a smooth mapping s : M — > V such 
that s(x) G V x for all The sections of V are to be 
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multiplied in the exterior sense, as they represent anti- 
commuting degrees of freedom; hence the proper object 
to consider is the exterior bundle, AV. 

It is a beautiful fact that there exists a unique ac- 
tion of the Lie superalgebra g on the sections of AV 
by first-order differential operators, or derivations for 
short. (Be advised however that this canonical g-action 
is not well-known in physics or mathematics.) 

The manifold M is a symmetric space, i.e. a Rieman- 
nian manifold with G-invariant geometry. Its metric 
tensor, g, uniquely extends to a second-rank tensor field 
(still denoted by g) which maps pairs of derivations of 
AV to sections of AV, and is invariant with respect to 
the g-action. This collection of objects — the symmet- 
ric space M, the exterior bundle AV over it, the action 
of the Lie superalgebra g on the sections of AV, and 
the g-invariant second-rank tensor g — form what the 
author calls a Riemannian symmetric superspace, ^ . 

7.3. Non-linear sigma model. — According to the 
Landau-Ginzburg-Wilson paradigm of the theory of pha- 
se transitions, the large-scale physics of a statistical me- 
chanical system near criticality is expected to be con- 
trolled by an effective field theory for the long-wave 
length excitations of the order parameter of the system. 

Wegner is credited for the profound insight that the 
LGW paradigm applies to the random-matrix situation 
at hand, with the role of the order parameter being taken 
by the matrix Q. He argued that transport observables 
(such as the electrical conductivity) are governed by 
slow spatial variations of the <2-field inside the saddle- 
point manifold. Efetov skilfully implemented this in- 
sight in a supersymmetric variant of Wegner's method. 

While the direct construction of the effective contin- 
uum field theory by gradient expansion of (flU l is not 
an entirely easy task, the outcome of the calculation 
is pre-determined by symmetry. On general grounds, 
the effective field theory has to be a non-linear sigma 
model for the Goldstone bosons and fermions of if 
{0 A } are local coordinates for the bundle V with metric 
gAB(ty) = g(d / dty A ,d / dty B ), the action functional is 

S = aJ d d xd^ A g AB ($)d^ B . 

The coupling parameter a has the physical meaning of 
bare (i.e. unrenormalized) conductivity. In the present 
model a = NW 2 a 2 ~ d , where W is essentially the width 
of the band random matrix in units of the lattice spacing 
a (the short-distance cutoff of the continuum field the- 
ory). S is the effective action in the limit z\ = Zi- F° r a 
finite frequency co = zi — Zi, a symmetry-breaking term 
of the form icov f d d xf((j)), where v = N(7tX)~ l a~ d is 
the local density of states, has to be added to S. 



By perturbative renormalization group analysis, i.e. 
by integrating out the rapid field fluctuations, one finds 
for d = 2 that a decreases on increasing the cutoff a. 
This property is referred to as asymptotic freedom in 
field theory. On its basis one expects exponentially de- 
caying correlations, and hence localization of all states, 
in two dimensions. However, a mathematical proof of 
this conjecture is not currently available. 

In three dimensions and for a sufficiently large bare 
conductivity, the renormalization flow goes toward the 
metallic fixed point (a — ► °°), where G-symmetry is 
broken spontaneously. A rigorous proof of this impor- 
tant conjecture (existence of disordered metals in three 
space dimensions) is not available either. 

7.4. Zero-mode approximation. — For a system in 
a box of linear size L, the cost of exciting fluctuations 
in the sigma model field is estimated as the Thouless 
energy Ejh = a/vL 2 . In the limit of small frequency, 
| CO | -C ZiTh, the physical behavior is dominated by the 
constant modes <^ A {x) = A (independent of x). By 
computing the integral over these modes, Efetov found 
the energy-level correlations in the small-frequency limit 
to be those of the Gaussian Unitary Ensemble. 
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